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We propose a scheme which can realize an extended two-component Bose-Hubbard model using 
polaritons confined in an array of optical cavities. In addition to the density-dependent interactions, 
this model also contains nonlinear coupling terms between the two components of the polariton. 
Using a mean-field calculation, we obtain the phase diagram which shows how these terms affect 
the transition between the Mott insulator and the superfiuid phase. In addition, we employ both 
a perturbation approach and an exact diagonalization method to gain more insights into the phase 
diagram. 
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I. INTRODUCTION 

Over the past decade, with the development of the 
technique of the optical lattices, one could manipulate 
cold atoms in the periodic potentials, which serve as 
an ideal platform to explore various many-body models 
in condensed matter physics [H-Q. As one of the im¬ 
portant fundamental models in many-body physics, the 
Bose-Hubbard model(BHM) has attracted lots of atten¬ 
tion since it was first introduced in 1963 Q. In the single 
component BHM, there exists a quantum phase transi¬ 
tion between the Mott insulator (MI) phase and the su¬ 
perfiuid (SF) phase [1, @]. This transition has been ex¬ 
plored in recent cold atom experiments Q . When this is 
extended to two components of bosonic modes, the phase 
diagram becomes more complicated correspondingly [3- 
1 - 

The standard single-component BHM only includes the 
on-site repulsive interaction U and the tunnelling t be¬ 
tween nearest neighbours. The interaction term favors 
the MI phase in which each site hosts a definite number 
of particles; whereas the tunnelling term favors the su¬ 
perfiuid phase where the particles are delocalized. The 
competition between these two terms drives a quantum 
phase transition between these two phases. For the two- 
component BHM, besides the intra-species interaction Ui 
and the tunnelling ti for each component (where * = 1,2 
labels the components), there is the additional inter¬ 
species interaction between the two components [l|. In 
recent years, a lot of attention has been paid to the ex¬ 
tended single (two or more) component BHM, and one 
kind of the extended BHM has additional tunnelling be¬ 
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tween the next-nearest neighbours or long range interac¬ 
tions 0, [H [H® , while another kind has additional 
nonlinear coupling between different components [iM3- 
Expectedly, these additional terms in the extended BHM 
give rise to a more complicated phase diagram pin- 


In a previous work we investigated a system of 
two-component polariton confined in a single optical mi¬ 
crocavity. Polariton is a kind of bosonic quasiparticle in 
the light-matter system, and Bose-Einstein condensation 
(BEG) of microcaviW polaritons have been realized in ex¬ 
periment recently . In this model, in addition to 

the density-density interaction (Kerr nonlinearity), there 
also exist two types of nonlinear coupling between the 
polariton components. In the current work, we extend 
this model by consider a one dimensional array of such a 
cavity polariton system, which realizes a two-component 
BHM. We will calculate the phase diagram of the sys¬ 
tem and focus on how the inter-species interactions and 
the nonlinear coupling terms will affect the transition be¬ 
tween the MI phase and the SF phase. 


This paper is organized as follows. In Sec. H we present 
the model Hamiltonian that describes an extended two- 
component Bose-Hubbard model. In Sec. HI, we provide 
the ground state phase diagram obtained from the mean- 
field decoupling approach. To gain futher insights and 
examine the effects of different terms in the Hamiltonian, 
we present a perturbative analysis in Sec. IV. In Sec. V, 
we give a beyond-mean-field exact diagonalization calcu¬ 
lation. In particular, we calculate the number fluctuation 
per cavity, and establish connections between this result 
with the mean-field result presented earlier. Finally, a 
summary is presented in Sec. VI. 
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FIG. 1: (Color online) (a) An illustration of a one-dimension 
array of optical microcavities. Each cavity host a system of 
two-component polaritons. (b) The atomic level structure and 
coupling laser fields that give rise to the two-component po- 
lariton system in each cavity. Here Qij is the coupling between 
the dipole transition (|i) •<->■ |j)) induced by the cavity mode, 
Q.ij and are the Rabi frequencies of the driving laser fields 
and driving microwave fields respectively. For details, please 
see Ref. 1^ . 


II. MODEL 


Our underlying system is schematically shown in 
Fig. [TJa). We consider a one-dimensional array of con¬ 
nected optical cavities, with the connection provided by 
the photon tunneling between neighboring cavities. In¬ 
sider each cavity, we have an ensemble of bosonic atoms 
whose level structure is sketched in Fig. [TJb). We la¬ 
bel the atomic hyperfine states as |i) with z = 1, 2,..., 7. 
Three of these states (states |1), |2), and |3)) belong to 
the electronic ground manifold, and the other four (states 
14), 15), 16), and |7)) belong to the electronic excited 
manifold. The ground states |1), |2), and |3) are dipole 
coupled to the excited states |4), |5), and |6), respec¬ 
tively, by the cavity field, with corresponding coupling 
strengths gn, 325 and Additionally, states |2) and 
13) are coupled to |4) by external laser fields with coupling 
strengths fl 24 and 1234 , respectively. Finally, within the 
excited manifold, states |5) and |6) are coupled to |7) by 
microwave fields with corresponding coupling strengths 
and fli/ 2 - e, <5, A, A 5 and Ag are various detunings 
between the driving field and the corresponding atomic 
transitions, as labelled in Fig. HKb). 

In the limit of weak excitation where the atomic popu¬ 
lation in the excited levels is negligible, we can construct 
two polariton modes for each cavity. The corresponding 
annihilaton operator for the two polariton modes in the 


zth cavity is given by 
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where ai is the cavity photon annihilation operator for 
the ith cavity, and 


1 

V ai 

S3. = _ 5 ]| 1 ),,( 3 |, 


are the collective atomic operator, with Nai being the 
total atom number, in the zth cavity. For simplicity, we 
have taken 1724 = ^^34 = gi = \/Nai gu, and Wi = 

y/gf + ^. The total Hamiltonian of the system reads: 
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+ T+[pUpIPu + PlP2i)P2. + pUpI.Pu + 4 ^ 2 ,)^!.] 

+ T-[pI{pIPu - pIP2.)P2. + pUpI.Pu - pIP2.)Pu]. 

( 2 ) 


represents the Hamiltonian of a two-component polariton 
in the zth cavity, the derivation of which can be found in 
our previous work (2lj |. and 

Hi.,,) =-t{PlPi, + PIP 2 ,} , (3) 

describes the tunneling of polaritons between adjcent 
cavities [ 2 ^. Here, 4(2)i ^ 1 ( 2 ). obey bosonic com¬ 
mutation relation: P^j] = 0, [Pj^p Pj^] = 0, and 

[Pci, Pp^] = Sa,i 3 Sij. As shown in Ref. [2l|, the key 
parameters in Eq. (ED, W, V 2 , U and can be tuned 
over a large extent by appropriately controlling the laser 
intensities and frequencies. Here to avoid instability, we 
only consider repulsive interactions such that Vi _2 and U 
are all positive. In addition, if we simultaneously change 
the signs of T^, the physics remains unchanged as that 
sign change can be absorbed by a redefinition of the po¬ 
lariton modes. Hence we will only consider the case with 
T'^ > 0, T~ < 0 for simplicity. 

In comparison to the single-component BHM, for 
which there are plenty of theoretical and experimental 
investigations, our model contains three new key param¬ 
eters: the inter-species on-site interaction characterized 
by the interaction strength U, the two nonlinear coupling 
terms characterized by the coupling strength T^. Our 
work will focus on elucidating the effects of these terms. 
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III. MEAN FIELD PHASE DIAGRAM 

Based on the mean field decoupling theory we 
introduce the superfluid order parameters: <^ 1 ( 2 ) = 

(Pi( 2 ).i). By the approximation 

= 4'1{2)Pi{ 2)J + <l^l{2)Pl(2),i ~ ?^ 1 ( 2 )^ 1 ( 2 ) > 

the total Hamiltonian can be decoupled into the following 
form: 

no = - ^^{p}Pl + p^P 2 ) + yPfA" + Y^^Pi 

+ UpIPiP^P2 

+ T+[PliPlPi + P^P2)P2 + pHpIPi + P2P2)Pi] 
+ t-[pI{pIpi - P^P2)P2 + pHpIpi - P2P2)Pi] 
+ + 02 (/> 2 ) , 
nt=- zt{(l)\Pi + + h.c.) , 

( 4 ) 

where /i is the chemical potential and z = 2 is the number 
of nearest neighbors. For simplicity, here, the site index 
is neglected. Since Hamiltonian 'Hq keeps local particle 
number conserved, in principle, we can numerically di¬ 
agonalize "Ho in a subspace with fixed total number of 
particles. 

As a reference, let us first set the nonlinear coupling 
strength = 0, and the tunneling rate t = 0. Un¬ 
der this situation T-Lq describes a two-component Bose 
gas with density-density interactions, and it preserves 
the number of polaritons in each component. The eigen¬ 
states therefore correspond to Fock states \ni,n 2 ) with 
definite integer values of ni and n 2 , and with correspond¬ 
ing ground state energy E(rii,n 2 ), where Ua {a = 1, 2) 
represents the number of polaritons in component a in 
each cavity. Note that due to the conservation of po- 
lariton numbers in individual component, we need to in¬ 
troduce two chemical potentials /ri and /r 2 for the two 
polariton modes. The ground state is determined by 
the relative size of Vi, V 2 and U. In Fig. [5] the ground 
state phase diagram is shown for the case with equal 
intra-species interaction strength, i.e., Vi = V 2 = V. If 
we focus on the situation with equal chemical potential 
(/ii = /r 2 = /r), it is quite straightforward to show that 
if the inter-species interaction is smaller than the intra¬ 
species one (i.e., U < V, see Fig. [UJa)), then the we 
have a single ground state with ni = n 2 for even total 
particle number, and two degenerate ground states with 
ni — 71-2 = ±I for odd total particle number. By contrast, 
iiU > V (see Fig. [2Kb)), no matter the particle number is 
even or odd, we have two degenerate ground states with 
ni = 0 or 712 = 0. 

The Fock state |7ii,7i2) discussed above corresponds 
to the Mott regime with 4>i = 4)2 = 0. As inter-cavity 
tunneling is turned on, (pa may take finite values and 
the system enters the superfluid regime. At zero tem¬ 
perature, the transition from the MI to the SF regimes 




FIG. 2: (Color online) The ground state of Ho with = 0 
and t = 0. Here we have to modify the first term in 'Ho to be 
—/iiP/Pi — fJ, 2 P 2 P 2 - The values in the bracket indicate those 
of ni and 712 , respectively. We have used Vi = V 2 = V, and 
the value of U is taken as U = 0.5 in (a), and U = 1.5 in (b). 
In all our figures, the quantities U, and are expressed 
in units of V. 

represents a quantum phase transition. Our goal is to 
map out a phase diagram by calculating the boundary 
between these two quantum phases. As we want to focus 
on the effects of the terms characterized by U and 
for most of our calculation we will choose Vi = V 2 = U 
as the units for energy, and consider the situation with 
U <V and U > V respectively while |T=*=| <C V. 

Figure O shows six examples of the phase diagram 
obtained from the mean-field decoupling approach with 
U <V in FigjSKa^c) and U > V in FigjSKd^f). We have 
also checked that essentially the same phase diagram can 
be obtained using the Gutzwiller method. The details of 
the phase diagram will be described in the following. 

A. Small Inter-species Interaction {U < V) 

In Fig.jSKa), the nonlinear coupling term is turned off, 
i.e., = 0. We again focus on the situation with equal 

chemical potential (/ii = /i 2 = /i). Here, the MI regime 
exhibits an “even-odd effect”, i.e., the MI region with odd 
occupation (the number of polaritons per cavity Titot = 
271-1-1) is smaller than its nearest neighbors with Titot = 
271 and 2(71 -|- 1). A qualitative interpretation of this 
effect can be provided by studying the excitation gap of 
a MI [ 1 ^. For even occupation MI with Titot = 27i, the 
excitation gap is given by 

^ 2 n = E{n, 71 -b 1) -V E{n — 1, 71 ) — 2E{n, n) = V, 

while that for Wot = 27i -|- 1 is given by 

A 2 „+i = E{n, n) + E{n -b 1, n + 1) — 2E{n, n + 1) = U . 

Since V > U, MI regime with even occupation will 
be more robust and hence exists in a larger parameter 
space. The excitation gap will compete with the tun¬ 
neling energy I]a=i. 2 (^a,i-P«b-i-i) which can be effec¬ 
tively strengthened if there exists degeneracy. As the 
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FIG. 3: (Color online) The mean-field phase diagram where 
the line represents the boundary between the MI and the SF 
regime. The region below the line is the MI regime, and the 
number indicate the total number of polaritons per cavity, 
ntot- The parameter used are (a) U = 0.5, T"*" = T~ = 0; 
(b) U = 0.5, T+ = 0, T" = -0.05; (c) U = 0.5, T+ = 0.05, 
T- = -0.05; (d) U = 1.5, T+ = T- = 0; (e) U = 1.5, 
T+ = 0, r- = -0.05; (f) U = 1.5, T+ = 0.05, T" = -0.05. 
All the quantities are in units of Vi = V 2 = V. 


filling number increases, the tunneling energy grows and 
the MI regimes with large occupation number become 
smaller and smaller. However, when the nonlinear cou¬ 
pling terms are present, i.e., ^ 0, this “even-odd 

effect” is weakened and may disappear completely as 
shown in Fig. Ec). An interpretation of effect is pro¬ 
vided by perturbative analysis in Sec. IV. Furthermore, 
the nonlinear coupling terms could change the degener¬ 
acy of the ground state which will be discussed in detail 
as follows. 

As mentioned previously, if t and are all zero, the 
ground states are doubly degenerate when ntot = ni -|-n 2 
is odd. If t is switched on, due to the super-exchange 
interaction induced by the inter-cavity tunnelling, the 
degeneracy would be lifted up and the ground state be¬ 
comes the superposition state of ni — n 2 = ±1 with 
equal weight [9|. This is indicated in Fig. 21 where we 
have defined the average particle numbers per cavity 
Ha = {PM {a = 1,2) and the fluctuations Aua = 
{{PlPc?) - Antot = {{pIPi + P^P2?) - 

(P/Pi -I- ^ 2 ^ 2 )^- For Fig. HDa) and (c), we take a rel¬ 
atively small chemical potential with n/V = 0.25, for 
which the average particle numbers are ni = n 2 = 0.5. 
In the absence of the tunneling, the ground state is dou¬ 
bly degenerate and are represented by the Fock states 
|1,0) and |0,1). For 0 < t < tc, where p is the criti¬ 
cal tunneling strength beyond which the system changes 
from MI to SF and is indicated in the figure by the verti¬ 


cal dashed line, the fluctuations of the total number per 
cavity Aritot vanishes, while Ani and An 2 are both hnite, 
as shown in Fig. |D)c). Therefore, we can conclude that 
the ground state is an equal-weight superposition state 
of 11,0) and |0,1). For the case depicted in Fig. UDb) and 
(d), we used a larger chemical potential fi/V = 1 and 
the average particle number is ntot = 2. In this case, 
the MI state is represented by the Fock state |1,1) with 
vanishing fluctuations in both Antot and Ani^ 2 - 

If T~ is turned on, the ground state for t < tc would 
be a superposition state Fock states |ni,n 2 ) with fixed 
ni + n 2 = ntot- This state is non-degenerate if ntot is 
even, and can be written as 

IG) = X! fni,n2\ni,n2 ), (5) 

ni,n2 


with = (-1)' ^2 fn 2 ,ni- An example is repre¬ 

sented in Fig. [SDa) and (c) with ntot = 2. For t < tc, 
we note that the particle number difference remains zero 
and ni ^2 have finite fluctuations while Antot vanishes, so 
the ground state is a superposition of states |1,1), |0, 2) 
and |2,0). By contrast, when ntot is odd (see Fig. [SDb) 
and (d)), T~ separates the whole Hilbert space into two 
degenerate subspaces with ni > n 2 and ni < n 2 , re¬ 
spectively, and the corresponding ground states have the 
form 


|G+) = c+,„Jni,n2), (6) 

ni >712 

and 

|G-) = c;:i.„Jni,n2), (7) 

ni <712 

TL 2 ^ ^ ^ ^ 1 

with = (-1) 2 One can readily ver¬ 

ify that the T~ term has vanishing matrix element be¬ 
tween |G+) and |G_). Thus the double degeneracy of the 
ground state is preserved even if T~ is turned on, and the 
“even-odd effect” remains, as shown in Fig. EKb). 

However, if we turn on T'^, the degeneracy will be 
lifted, as the T+ term could couple the two degenerate 
subspaces, i.e., it has non-vanishing matrix elements be¬ 
tween |G-|-) and |G_). This explains the weakening of 
the ‘“even-odd effect” as shown in Fig. [SDc). 


B. Large Inter-species Interaction {U > V) 

In Fig. HDb), we have known that the ground states, 
i.e., |0,'ntot) and |ntot,0), are doubly degenerate for both 
the even and the odd occupation when U > V, = 0 

and = p ,2 = M at f = 0. In Fig. [HI we plot the 
occupation numbers and their fluctuations as functions 
of tunneling rate t. we can see that, in the MI region 
(0 < t < tc), the population difference is ni — n 2 = ±n.tot 
while the number fluctuations Antot and Ani ^2 are all 
zero. It tells us that the ground state is either |0,ntot) 
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FIG. 4: (Color online) The average particle numbers and their 
fluctuations obtained by the Gutzwiller method. For all the 
subplots, we have Vi = V 2 = V, U = 0.5, and = 0. The 
chemical potential is taken as = 0.25 for (a) and (c); jJ. = 1 
for (b) and (d). The vertical dashed lines indicate the critical 
tunneling rate tc beyond which the system changes from MI 
to SF. 
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the tunneling energy in this nondegenerate state is much 
lower than that in |0, ntot) or |ntot, 0) while the excitation 
gap decrease is very small. Thus the MI region with even 
occupation is enlarged by T~ and the “even-odd effect” 
reappears. 

Finally, Fig.[3Kf) shows that the “even-odd effect” van¬ 
ishes again if T+ is also turned on. The explanation for 
this is quite similar to the case when U < V: the T+ 
term lifts the degeneracy for odd ntot as its matrix ele¬ 
ment between |C?+) and |(1?-) is non-zero. 

We comment in passing that, from Fig. [3Ka)~ (f), we 
may notice that the MI region with ntot = 1 remains 
unchanged for different values of T^. This is because 
the terms’ matrix elements are all zero in this space. 
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FIG. 6: (Color online) The average particle numbers and their 
fluctuations obtained by the Gutzwiller method. For all the 
subplots, we have Vi — V 2 = V, U = 1.5, and = 0. The 
chemical potential is taken as /r = 0.5 for (a) and (c); /t = 1.5 
for (b) and (d). The vertical dashed lines indicate the critical 
tunneling rate tc beyond which the system changes from MI 
to SF. 


FIG. 5: (Color online) The average particle numbers and their 
fluctuations obtained by the Gutzwiller method. For all the 
subplots, we have Vi = V 2 = V, U = 0.5, T”*" = 0, and 
T~ = —0.05. The chemical potential is taken as /r = 1 for 
(a) and (c); /r = 1.75 for (b) and (d). The vertical dashed 
lines indicate the critical tunneling rate tc beyond which the 
system changes from MI to SF. 


or |ntot 7 0), regardless whether ntot is even or odd. In 
addition, the excitation gap is independent of ntot, i-e., 
A = V. As a result, there is no “even-odd effect” as can 
be seen from Fig.[3Kd). 

However, when T~ is turned on, the “even-odd effect” 
emerges as can be seen in Fig. [SKe). This is because for 
finite T~, the ground state in the Mott insulator region 
is non-degenerate when ntot is even, and takes a similar 
form as in Eq. An example is shown in Fig. [TKa) 
and (c). By contrast, for odd ntot, T~ can preserve the 
ground state’s double degeneracy and the two degenerate 
ground states take similar forms as in Eqs. m and 0, 
but with It is easy to find that 


It is not difficult to notice from FigITKa) and (c) that 
there is a discontinuous jump in both the particle num¬ 
bers and their fluctuations at the critical tunneling rate 
tc- The jump in Antot, in particular, indicates that the 
phase transition from MI to SF may be of first order for 
even ntot when U > V and T~ ^ 0. To verify this, we 
fix UjV = 1.5, T+ = 0, T- jV = -0.05, /r/F = 1.5, 
and plot the ground state energy E{(j)i,(j) 2 ) with differ¬ 
ent t in Fig. [SI When t is very small, (j)i = (j )2 = 0 (MI) 
is the global minimum point (see Fig. ISKa)). As t in¬ 
creases, additional local minimum points with finite 4 )i ^2 
(metastable SF) arise, as shown in Fig. jS^b). If we con¬ 
tinue to increase t, Fig. He) displays that the MI state 
with (()i 2 = 0 becomes a local minimum point (i.e., a 
metastable MI). Finally, the local minimum point be¬ 
comes a local maximum point if t is sufficiently large as 
shown in Fig. |H{d). In conclusion, metastable SF and 
metastable MI states exist near the boundary of the MI 
lobes in Fig. |3](e) (between the dashed and solid line) 
while ntot is even. This existence of such metastable 
states is a tell-tale signature that the MI-SF phase tran- 
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FIG. 7: (Color online) The average particle numbers and their 
fluctuations obtained by the Gutzwiller method. For all the 
subplots, we have Vi = V2 = V, U = 1.5, T"*" = 0, and 
T~ = —0.05. The chemical potential is taken as — 1.5 for 
(a) and (c); /r = 2.5 for (b) and (d). The vertical dashed 
lines indicate the critical tunneling rate tc beyond which the 
system changes from MI to SF. 


sition in this case is of first-order. A similar situation can 
be found in spin-1 bosons which can also host the first- 
order phase transition [mil, but they are induced by 
different effects. In Fig. |H1 another feature is also worth 
our attention. Only one of the two species is dominant 
in the SF state, i.e. \4>i\ ^ \(j) 2 \ or \(j) 2 \ ^ \4>i\ and they 
are degenerate. This is in consistent with Fig. [Tic) and 
(d) which indicate that only one component’s fluctua¬ 
tion is prominent while the other one’s is close to zero. 
In comparison, both components have significant occu¬ 
pation fluctuations in SF state when U < 17, as shown in 
Fig. Sic) and (d), and Fig. jSlc) and (d). 

Furthermore, if T~ = 0 but yf 0, not only the 
even occupation’s but also the odd occupation’s MI-SF 
phase transition will be of first-order, as shown in Fig. [HI 
Because all the MI ground states’ degeneracy is lifted in 
the presence of T'^ regardless whether ritot is even or odd. 



FIG. 8: (Golor online) The energy contours obtained by di¬ 
agonalizing Eq. (|4]) for U = 1.5, /r = 1.5, = 0, and 

T~ = —0.05 (all in units of V\ = V2 = V). The tunneling rate 
is given by (a) t/V = 0.05; (b) t/V — 0.6; (c) t/V = 0.08; 
(d) t/v = 0.1. Though the order parameters, could 

also be complex numbers, their relative phase must be 0 or 
TT in the ground state of Eq. 0. Hence we can perform the 
calculations assuming (f )\^2 are real, without loss of generality. 


a 


< 



1.2 _-^^^^ 

1U»*0*0000 

. An, 

0.4 * 

02V^ntot 

0.0 TVWVVW" 

0.00 0.05 0.10 0.15 0.20 

t/v 


jc) 

>• 0*0000 

• Ani 
■ An, 

V An,ot 

TVVVWW" ■■■■••■•■ ■•• 








3.0 

2.5 
2.0 

1.5 
1.0 
0.5 
0.0 

0.00 0.05 0.10 0.15 0.20 

2.5 

2.0 

1.5 
1.0 
0.5 
0.0 

0.00 0.05 0.10 0.15 0.20 

t/v 


• Ani 

■ An2 
V Antot' 


(b) ' ♦•♦*o*o*o*oco«q 

O «1 ♦ «2 




IV. PERTURBATIVE ANALYSIS 

To gain more physical insights into the phase dia¬ 
gram, we consider the inter-component interaction (the 
U-term) and the nonlinear coupling (the T^-terms) as 
perturbations and examine how they affect the phase di¬ 
agram. Furthermore, in order to distinguish the contri¬ 
bution from each of them, we have shown how each of 
these nonlinear perturbations affect the mean-field phase 
diagrams in Fig. [101 

In Fig. [Tol a), we show the “unperturbed” phase di¬ 
agram by taking U = 0 and = 0. Here we sim¬ 
ply have two uncoupled single-component Bose-Hubbard 
model. Since we have chosen the intra-component inter¬ 
action strength VT = V 2 = the two individual phase 
diagrams completely overlap with each other. In partic- 


FIG. 9: (Golor online) The average particle numbers and their 
fluctuations obtained by the Gutzwiller method. For all the 
subplots, we have Vi = V2 = V, U = 1.5, T'^ = 0.05, and 
T~ = 0. The chemical potential is taken as = 1.5 for (a) 
and (c); fi = 2.5 for (b) and (d). 


ular, each MI region is characterized by equal number of 
the two polariton modes, and hence only even TVtot MI 
regions are present. Along the horizontal axis {t = 0), 
the MI region with Ntot = 2n occupies a region between 
{n — 1)17 < ^ < nl7. At the boundary ^ = nV, the four 
Mott states |n, n), |n, n -I-1), |n -|- 1, n) and |n -|- 1, n -I- I) 
are energetically degenerate with energy —n{n -|- 1)17. 

In Fig. nUT b). we show the phase diagram in the pres¬ 
ence of a small inter-component intereaction f/ <C 17. As 
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FIG. 10: (Color online) The mean field phase diagrams for 
different nonlinear pertuibative terms. Here, Vi = V 2 = V. 
(a) phase diagram with U = 0 and = 0; (b) phase diagram 
with U = y/30 and = 0; (c) phase diagram with = 
H/eO and U = T- =0. 


one can see, the [/-term induces MI regions with odd 
n-tot = 2n -b 1, which for small [/, occurs near the bound¬ 
aries between two even MI regions at n = nV. At t = 0, 
these odd MI regions are represented by two degenerate 
Fock states \n,n + 1) and |n -I- l,n), and occupy a re¬ 
gion with /i £ The values of /i< and /r> can 

be readily obtained from i/(n,n,/r<) = E{n,n + l,/r<) 
and E{n, n + 1, /r>) = E{n -I- I, n -I- I, ir>), from which 
we obtain the width of the odd MI region in the /i-axis 
as — /i< = U, which is exactly the excitation gap 
^ 2 n+i of the MI state. The above analysis provides a 
more quantitative argument for the even-odd effect men¬ 
tioned earlier. Furthermore, since the tunneling energy is 
roughly proportional to a/ n{n + 1), the MI regimes with 
occupation number 2n -I- 1 become smaller as n grows. 
Our mean-field numerical results, as shown in the inset 
of Fig. [TUf b). are in complete agreement with this anal¬ 
ysis. 

In Fig. ITOl cl. we show the phase diagram in the pres¬ 
ence of the nonlinear coupling T+-term, while keeping 
U = 0 and T~ = 0. In this case, in the absence of tunnel¬ 
ing, the conservation of the number of polaritons in each 
component is broken, while the total polariton number 
ritot for each cavity remains conserved for Hamiltonian 


EIq. We also observe the appearance of odd MI regions 
except for ritot = 1- Furthermore, the odd MI regions 
grow in size as ritot increases, which is in stark contrast 
with the previous case. These properties can be intu¬ 
itively understood as follows. Consider an odd MI region 
with ritot = 2n + 1 at t = 0, which is characterized by 
two degenerate Fock states |n,n + 1) and |n + l,n) in 
the absence of T+. When a small T+ is turned on, a 
direct coupling between these two states is induced with 
the corresponding matrix element given by 2n(n + 

To first order in T+, the ground state energy is shifted 
down by this amount. The presence of T+-term will also 
lower the energies of even MI states, but the energy shift 
is quadractic in T+ as the unperturbed even MI states 
are nondegenerate. This explains the appearance of odd 
MI regions with ritot = 2n -|- 1 for n > 0. There is no MI 
region for Wot = 1, because transition matrix element 
between Fock states |0,1) and |1,0) vanishes. In addi¬ 
tion, it is not difficult to show that the width of the odd 
MI region on the /i-axis is 0(4n(n + I)T+), which grows 
when n increases, in good agreement with the numerical 
results. 

Actually, no matter U is zero or not, we can always 
use the above picture, i.e., first order energy shift to the 
state with degeneracy and second order energy shift to 
the state without degeneracy, to understand the effect 
of T+. And it also applies to the situation even U > 
V . As mentioned previously, the ground states with odd 
occupations are doubly degenerate when U > V, T~ 7/ 
0 at t = 0, and T~^ has nonvanishing matrix elements 
between them which can induce a first order energy shift. 
At the same time, the energy shift is second order in the 
even occupation case. On the other hand, the tunneling 
energy grows as T+ increases. Therefore, the “even-odd 
effect” disappears in Fig. |3Kf). 

Finally, let us consider the effect of the T“-term. In 
the case U = 0 and r+ = 0, the T^-term alone would 
not induce odd MI regions. This is due to the fact that, 
unlike the T+-term, this term does not induce a direct 
coupling between |n, n + 1) and |n -|- 1, n) states, i.e., the 
corresponding matrix element vanishes. Hence a small 
T~ does not have any noticeable effects on the phase di¬ 
agram. However, the T^-term can be regarded as a non¬ 
linear coupling between the two components of the po¬ 
lariton, whose effective sign depends on the population 
difference. As a result, in the superfluid region, when 
the relative population rii — 112 changes sign, the rela¬ 
tive phase between the two order parameters (j)i and (j )2 
will change from 0 to tt. As this effect is already present 
in the single-cavity system which we investigated in de¬ 
tail in Ref. [2l|, we do not provide a detailed discussion 
here. Instead, we just present the phase diagram for the 
superfluid region in Fig. 1111 in which OSF and ttSF are 
superfluid phases with relative phase between tpi and 4>2 
being 0 and tt, respectively. The transition between OSF 
and ttSF is of first order. 

The perturbation calculation presented above allows us 
to better understand the phase diagram shown in Fig. |31 
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FIG. 11: (Color online) The phase diagram for the first order 
phase transition shown in the 3D parameter space (, 

The blue surface is the phase boundary, while 
the OSF and ttSF is below and above it respectively. 

In particular, we can now explain why the inter-species 
interaction give rise to the even-odd effect, and how the 
nonlinear coupling term T~^ weakens the even-odd effect. 


V. EXACT DIAGONALIZATION 

So far we have investigated a homogeneous system of 
inter-connected cavity in ID using a mean-field approach. 
In this part, the exact diagonalization method is used to 
study this model. To make the calculation manageable, 
we consider a finite number, Nc, cavities with Np to¬ 
tal particles. The whole Hilbert space are spanned by 
the Fock state basis {ipi) = - where 

denotes the number of polaritons in component-a 
in fcth cavity, and they are constrained as 0 and 

= ^p- The dimension of the Hilbert space 
is therefore Dm = write the Hamiltonian 

into a matrix form under this Fock state basis using the 
periodic boundary condition, and obtain the ground state 
(IG)) of this large sparse matrix through exact diagonal¬ 
ization. In Fig. 1121 we plot the total number fluctuation 
per cavity which is defined as 

Ari = (G|(n'=)2|G)-(G|ri'=|G)2 

where = n\ + is the total number operator for 
kth cavity. Under the periodic boundary condition, this 
quantity is indepedent of the cavity index k. We vary Np 
and Nc to some extent while restricting their ratio Np /Nc 
(i.e., the number of polaritons per cavity) to be 1, 2, or 3. 
As can be seen, the behavior of the number fluctuation is 
sensitive to Np/Nc- The parameters of Fig. fTW a) are the 
same as those used in Fig. fTHT c). Here one can see that 
in the limit t —>■ 0, the number fluctuation for systems 
with Np/Nc = 1 remains finite, which indicates the lack 
of MI region for Wot = 1- By contrast, in the same limit, 
the number fluctuations for Np /Nc = 2 and 3 vanish and 


more specifically An for Np/Nc = 3 tends to zero with 
a much steeper slope. We thus expect to see a large MI 
region for Wot = 2 and a small MI region for Wot = 3. 
All these are fully consistent with the mean-field phase 
diagram shown in Fig. w. The parameters used for 
Fig. fT^ b) are the same as those used for Fig. fTUT b). Here 
An for all three ratios or Np/Nc vanish in the limit t —>■ 0. 
The slopes of An near this limit indicate that two small 
MI regions for Wot = 1 and 3, and a large MI region 
for Wot = 2 in the mean-field limit is expected, which is 
again consistent with the results obtained earlier. 



FIG. 12: (Golor online) The fluctuation of the total occu¬ 
pation per cavity An, obtained from the exact diagonaliza¬ 
tion method. Here Vi = V 2 = V and in (a) U = T~ = 0, 
r+ = U/60; in (b) U = U/30, T+ = T~ = 0. The parameters 
used in (a) and (b) are the same as those used in Fig. llOl cl 
and (b), respectively. 


VI. SUMMARY 

In summary, we have presented a scheme to realize 
a two-component BHM with nonlinear inter-component 
coupling in a system of cavity polaritons. We mapped 
out the phase diagram showing the boundaries between 
the MI phase and the SF phase. Using several differ¬ 
ent approaches — the mean-field decoupling method, 
the Gutzwiller method, the perturbation calculation and 
the exact diagonalization, we show how the inter-species 
interaction and the nonlinear coupling terms affect the 
phase diagram, and particularly how they induce the 
first-order MI-SF phase transition and give rise to or 
weaken the “even-odd effect”. Additionally, the competi¬ 
tion between the nonlinear coupling strengths and T~ 
can drive a first order quantum phase transition within 
the SF regime that changes the relative condensate phase 
of the two polariton components. Through our study, 
we have obtained a clear understanding about this two- 
component BHM. In the future, we could also realize 
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an two-component BHM that breaks the time reversal 
symmetry [l6l Il7j| by manipulating the external fields’ 
relative phase. This model may host more exotic p-wave 
superfluid phases [l5l| . 
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